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SUMMARY 

An Improved finite-difference procedure based 
upon second order central differences Is developed 
In this paper. Several difficulties encountered 
in earlier works with fictitious stations that 
arise in using second order central differences, 
are eliminated by developing certain recursive 
relations. The need for forward or backward dif- 
ferences at the beam boundaries or other similar 
procedures is eliminated in the present theory. 

By using this improved theory, the vibration char- 
acteristics of pretwisted and tapered blades are 
calculated. Results of the second order theory 
are compared with published theoretical and exper- 
imental results and are found to be in good agree- 
ment. The present method generally produces close 
lower bound solutions and shows fast convergence. 
Thus, extrapolation procedures that are customary 
with first order finite-difference methods are 
unnecessary. Furthermore, the computational time 
and effort needed for this improved method are 
almost the same as required for the conventional 
first order finite-difference approach. 

INTRODUCTION 

Determination of the vibration characteristics 
of turbomachine blading by theoretical means has 
become very important in recent years. Blade 
failures are still not uncommon even with the 
present d'y material technology and improved design 
methodology. These failures are normally attribu- 
ted to fatigue which occurs when blades vibrate at 
or nerr resonant conditions. It thus becomes 
imperative that the designer be able to determine 
the natural frequencies and mode shapes of the 
vibrating blades as accurately as possible at the 
early stages of design. Several methods of solu- 
tion have been developed in the past in which only 
a few of the parameters such as pretwist, taper, 
asymmetry of cross section, disc and root flexi- 
bility, or of the construction complexities such 
as lacing wires and shrouding, are taken into 
account together. Inclusion of shear deflection, 
rotary inertia, warping, thermal and Coriolis 
effects in the vibration analysis becomes important 
for short and stubby beams subjected to high 
centrifugal force field. Recent research in heli- 
copter rotor blade and advanced turboprop blade 
vibrations reveals the necessity of including tne 
geometric nonlinearities up to varying degrees in 
the analysis in order to obtain fair prediction of 
the vibration characteristics and flutter boundar- 
ies. Some recent contributions in deriving the 
equations of motion include those of Kaza and 
Kvaternik [1], Hodges and Dowell [2] who used 
geometric-nonl inear theory and that of 
Subrahmanyam, et.al. [3] who used a linear theory 
but accounted for shear and rotary inertia effects, 
suitable for turbomachine blades. 
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The various theoretical methods of solution 
can be broadly classified as belonging to either 
the continuum model approach or the discrete model 
approach. In the continuum model approach, appli- 
cations of the potential energy, complementary 
energy or the Galerkin methods are well known. 

These methods produce upper bound solutions. The 
Reissner and the Dean and Plass methods which are 
classified as mixed methods [4, 5] were shown to 
produce upper bound solutions with relatively 
quicker convergence as compared to the classical 
approaches. In certain instances, the Reissner 
method may produce oscillatory convergence, depend- 
ing on the choice of shape functions used [6]. 

While the approximate methods mentioned above have 
certain advantages, at least one disadvantage is 
associated with the numerical evaluation of the 
integrals arising in the formulation of the fre- 
qi - ' equation. 

In the discrete model approach, the Holzer- 
Myklestad, Stodola, polynomial frequency equation, 
transformation. Station function and finite- 
element methods are well developed. Application 
of the Galerkin, finite-difference and collection 
method' to solve the equations of motion has been 
knc.' or quite some time. In general, the dis- 
crete model approaches produce lower-bound solu- 
tions due to the discretization of the distributed 
mass and elasticity. In certain finite-element 
applications, the convergence may be upper-bound 
depending upon the mass matrix formulation. Among 
the methods producing lower bound solutions, the 
first order finite-difference method is perhaps 
the single method that has attracted the greatest 
attention. Almost all the works dealing with the 
finite-difference method point out that the con- 
ventional first-order finite-difference method 
converges relatively slowly with mesh refinement. 
The improvement in accuracy expected with refine- 
ment of mesh size may be completely overshadowed 
by round-off and truncation errors that result 
from the increased matrix sizes. To avoid these 
difficulties, the Richardson’s extrapolation pro- 
cedure based on two or three solutions with dif- 
ferent mesh sizes [7, 8] has been applied. The 
extrapolated result shows improved accuracy but 
the results being extrapolated must possess a 
monotonic convergence. Furthermore, the extra- 
polated result may not necessarily be a bound. 

Relatively few works exist which use higher- 
order finite-difference techniques. Greenwood [9] 
used first-order and second-order finite differ- 
ences to analyze uniform cantilever beams in flex- 
ure, and one case with uniform breadth but varying 
depth. The fourth order differential equation was 
transformed into four first order equations, and 
the slopes and shearing forces were evaluated at 
half integer stations, while the deflections and 
bending moments were evaluated at integer stations. 
A central difference approximation of second order 
was used. One sided approximations were used to 
satisfy the boundary conditions. As an alternative 
approach, a complicated, but symmetric, method of 
applying the boundary conditions was illustrated. 
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It was concluded that one sided approximation gave 
more accurate results. It Is Interesting to note 
from Greenwood's results that for a tapered beam, 
the first-order central -difference theory produces 
more accurate results than the second-order theory. 
The reason was attributed to the lumping of non- 
uniform mass In the second-order theory. 


second moments of area about the axes specified by 
the subscripts. The breadth and depth taper param- 
eters are given by 

8 * < a 0 - b o ,/b o* 4 ■ < d o - c o )/c o < 9 ’ 10 ) 
Furthermore, 


In this paper a second-order central differ- 
ence approach is presented which eliminates most 
of the shortcomings discussed above. Uniform and 
tapered beams with or without pretwist are analyzed 
by using both the firstand second-order finite 
differences. The fictitious stations encountered 
in the development of the second-order theory are 
eliminated by using certain recursive relations 
derived by extending the central difference expres- 
sions for the boundary conditions given by the 
first-order theory. This logical extension In 
developing the recursive relations is shown to 
produce accurate results which converge more 
rapidly than results from the first-order theory. 
Furthermore, close lower-bound solutions can be 
obtained by the present approach, and extrapola- 
tions, such as are customary with the first-order 
theory, are not necessary when uncoupled vibra- 
tions are considered. 


n - z/L, ^ (11,12) 

For an untwisted blade, v - 0 and thus equations 
(1) and (2) are uncoupled. 

By making use of equations (3) to (12) in 
equations (1) and (2), performing the differentia- 
tions successively and grouping terms, the equa- 
tions of ,iotion can be written in the following 
form. 



+ f x pjy - 0 (13) 

dn 


VIBRATION ANALYSIS OF A PRETWISTEO 

TAPWM BLADE 


Equations of Motion 


Figure 1 shows a uniformly pretwisted and 
tapered blade of rectangular cross-section and the 
coordinate axes. The equations of motion [3, 10] 
for the coupled bending-bending vibrations of such 
a cantilever beam can be si own to be of the follow- 
ing form. 



d‘ 


dz 



* $}- ■ 0 


(1) 

(?) 


The sectional properties for a rectangular cross- 
section blade are 


! xx - I yy sin ? vn + i xx cos 2 vn. 


2 2 

I yy « Iyy cos vn + sin vn 


*YY ~ ! XX 


xy c 

I Y y- I oy (1 + «n)(l + Bn 


sin (2vn) 
3 


(3.4) 


(5.6) 


dn dn dn^ dn^ 


dn 5 


d 2 x 2 

+ v-T-^Y 

dn 


(14) 


In the above equations, \ - D A 0 L 4 /EI 0X and 

the coefficients a, b,...,f, q, r v are all 

functions of n* For brevity these functions are 
not presented here but can be found in reference 
[ 11 ]. 


Finite-Difference Equations for Derivatives 

First-Order Central Differences - The first 
order central differences for the derivatives of a 
function a at any arbitrary station i are 
given by 

♦i * k f*i-l + ♦ 1 +1 } H5) 

♦i ""7 {*1-1 ’*1 + *i+l} ( lfi ) 

♦l “ ^5 £» 1-2 + * i-i 1+1 + ♦ ,>2} 

♦i V *^r{*i-2 ~ S-l + (is) 


Ixx " *ox (1 + 4n) 3 (l + Bn), 
A - A q (1 "* Sn)(l + Bn) 


Second-Order Central Differences - The second 
order central differences for the derivatives of a 
. , function + at any arbitrary station i are 

('.81 given by 


In equations (1) to (8), E is the Youngs modulus, 
p the mass density of the blade material, p n 
the natural radian frequency, x and y the 
dynamic displacements of the centroid in xzand 
yz-planes, A the area at any section, A 0 the 
area of the root section, v the angle of pre- 
twist at the blade tip and I xx , I oy , etc., the 


♦i “ m{*i-2 ~ + 8»i + i -♦i^J (W) 

42 f * 1 " 2 + 1 ^ i - 1 ~ 3( ^' + 16>i+I 

(20) 
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♦i “ - ®»1-2 + 13 »1-1 ~ 13 »1+1 

+ *i+2~*1+ 3 j (21J 

♦N^{“*i- 3 + 12 »i-2- 3 *i-l + S6 »1 

-39 > . + i + 1 a >i+ 2-4i + 3} (22) 

in equations (15) to (22), the subscripts i-3, 

1-2, 1-1, 1, 1+1, 1+2 and 1+3 represent stations 
separated by the Interval h , where h - L/n, n 
being the number of segments Into which the beam 
of length L is divided. 

8oundary Conditions 

Tne boundary conditions associated with the 
coupled bending-bending vibrations of a pretwisted 
cantilever beam fixed at n » 0 and free at n ■ 1 
reduce to 

i I 

y « x «• y m x • 0 at d < o; 

II II III III 

y « x « y -x . 0 at i) ■ 1 (23,24) 


By using the central difference relations for the 
derivatives of x and y in the forms shown by equa- 
tions (15) to (13), one can easily show that 


y -l * y y X -1 - x l’ x o • >0 “ 0 
Vl - 2y n - Vl- x n+l - 2x n ' Vl 


y n+2 = 4y n - 4y . 


n-1 
x 


r n-2’ 


h+2 


4 x - 4x . + x _ 
n n-1 n-2 


(25) 

(26) 


(27) 


Equations (25) to (27) eliminate all the fictitious 
stations encountered in the development of the 
finite-difference formulation with first order 
central differences. In order to eliminate the 
fictitious stations y_? and y n +3 encountered 
in the second order central differences, one can 
assume a symmetry condition around the built-in 
end so that 

y_ 2 * 1^2 * * Xp ( 28 ) 


This assumption satisfies the boundary conditions 
(23) in terms of second order central differences. 
Equations (26) to (27) are rewritten in the fol- 
lowing form to establish a possible recursive 
relation: 


y n+l ~ y n-l = 2 ^ y n _ y n-l^’ 

Vl - x n-l * 2 < x n - x n-l> (29) 

y n+2 ' y n-2 “ 4 ^ y n ” y n-l); 

x n+2 " x n-2 ' 4(x n ' x n-l ) (30) 

and, by recursion 

y n+3 - y n-3 ' 6(y n " Vl ): 

X n+3 " x n-3 ” 6(x n " x n-l * < 31) 


Equations (31) can be used to eliminate y^ and 
x n+3* The error introduced by these assumptions 

can be evaluated by introducing the values of 
y n+l’ y n+2* and y n+3 obtained froni equations (29) 
to (31) into the second order finite difference 

II III 

expressions for y n and y n . It can be seen 

III II 

that y n « 0 •-.■Itile the expression for y n 
yields 

(»> 

Equation (32) states that the deflection at the 
(n-l) th station is the average of the deflections 
at the preceding and the succeeding stations. 

Thus, the deflection curve near the tip of the 
cantilever beam assumes a straight line form. 

Since the bending moment at the free end must be 
zero, the condition of constant slope near the tip 
is justified and, thus, the boundary conditions 
represented by equations (28) to (32) should give 
accurate results for a suitably large value of n. 

Finite-Difference Implementation 

8y introducing the finite-difference expres- 
sions for the derivatives, as given by equations 
(15) to (18) or (19) to (22) in terms of the func- 
tions x and y, into the differential equations 
(13) and (14), one obtains a coupled set of equa- 
tions for any arbitrary station. Each of this set 
can be evaluated for each station with i replaced 
by 1, 2....n successively. A system of 2n equa- 
tions can be written in terms of the variables 

yi. y2*--. yn. *1. ”2 *n- The 

parameters y_j, y_?, y n +i, y n +?, y+3, 

..... etc., can be eliminated By using equations 
(25) to (27) or (28) to (32). The resulting equa- 
tions can be represented in the familiar form of 
the eigenvalue problem 



" p n x 


(33) 


For the special case of zero pretwist, equation 
(33) reduces to two independent equations repre- 
senting uncoupled vibrations in the two principal 
planes as follows: 

[fi - P 2 y i] {yj- 0, [e - p 2 xj {x.J« 0 (34,35) 

where p and p are the uncoupled natural 
x y 

radian frequencies in the stiff xz-plane and flex- 
ible yz-plane respectively. 

In the preceding equations, £, £ and are 
square matrices of order (n x n) each, and { y } • 

and (x-j) are column matrices containing the 
displacement vectors of the n-stations. The eigen- 
values and the associated eigenvectors can be 
determined by using standard solution procedures. 

It may be noted here that the first and second 
order central difference methods generate frequency 
determinants of the same size for a given number of 
assumed stations. Both methods produce matrices 
that are banded. Each submatrix fl, S, C or Q 

possesses a band width of five for the first-order 


3 


i 
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theory and seven for second-order theory. Further- 
more, the matrices are non-symmetric for both 
methods developed. To save space the elements of 
the matrices are not presented here. 

NUMERICAL RESULTS AND DISCUSSION 

A FORTRAN computer program was developed to 
solve the eigenvalue problems defined by equations 
(31) to (35). The program was run on an IBM 370 
computer at the NASA Lewis Research Center. Use 
was made of an IMSL routine EIG2F [12]. This li- 
brary subroutine evaluates all the natural frequen- 
cies and the associated mode shapes (eigenvectors). 
As expected for a conservative system, it was found 
that all the eigenvalues and the corresponding 
eigenvectors are real for both the first and second 
order finite-difference applications. Several 
configurations of uniform and pretwisted blades 
were solved with and without taper. The results 
are presented below. 

Untwisted Uniform and Tapered Blades 

The following numerical data [13 to 14] were 
used to study the uncoupled vibrations of canti- 
lever beams of length 0.254 m with various breadth 
md depth taper ratios. The thickness ratio at 
the root section was kept as unity (square cross 
section), and the beams had an aspect ratio (L/b 0 ) 
of 40. A convergence study was made for a uniform 
beam divided into odd and even number segments. 

The frequency ratios presented in Table 1 indicate 
that the convergence is monotonic from below for 
both the first and second order central difference 
methods and that the convergence is continuous for 
even or odd values of n. Furthermore, the con- 
vergence shown by the second order theory is very 
rapid, and the lowest five flexural frequencies, 
though only the first three modes are shown in 
Table 1, can be obtained to within 0.195 percent 
error with n = 30. The first order theory shews 
errors of the order of 2.6 percent with n * 30, 
but an extrapolation procedure could be used to 
improve the accuracy. 

The natural frequencies and mode shapes are 
also calculated for tapered cantilever beams with 
various combinations of Lreadth and depth taper 
parameters. The results obtained by using the 
second order theory with n « 30 are in good agree- 
ment with the theoretical and experimental results 
presented in References [11, 13]. Because of space 
limitations, these results together with the re- 
sults for uncoupled axial and torsional vibrations 
of tapered cantilever beams were presented in 
Reference [14], 

Pretwisted Uniform Blades 

Recent investigations of the vibrational 
characteristics of pretwisted uniform blades in- 
clude those of Slyper [15] using the Stodola 
method, Dawson [16] using the Rayleigh-Ritz method, 
Dawson and Carnegie [17] using a transformation 
technique, Dokumaci, Thomas and Carnegie [18] using 
matrix displacement analysis and Subrahmanyam, 
Kulkarni and Rao [4, 5] using the potential energy 
and Reissner methods. Numerical data for the 
present study were chosen from these references, 
and the coupled bending-bending vibration frequen- 
cies and mode shapes were determined by both the 
first and second order finite-difference methods. 


Convergence trends remained similar to those pre- 
sented earlier though the convergence rates were 
not as fast as observed for the uncoupled case. 
Using n - 30 in the second order theory, the fre- 
quency parameter ratios, (x/x 0 ), [where x 0 is 
the fundamental frequency parameter for a blade 
with zero taper (exact value)], for the coupled 
bending-bending vibrations of uniform cantilever 
beams were obtained for typical values of breadth 
to depth ratios and pretwists. A typical set of 
such results are presented in Table 2, where a 
comparison of the present results is made with 
extrapolated results obtained from the matrix dis- 
placement method [18] and the first order finite- 
difference method. It can be seen from these 
results that the second order finite-difference 
method produces fairly accurate natural frequencies 
without extrapolations. It is worth mentioning 
that an extrapolated result may have an improved 
accuracy but can be either above or below the exact 
value. The present results obtained by the second 
order central differences show a maximum variation 
of the order of about *0.3 percent from the extra- 
polated results. In order to see whether the mode 
shapes for the coupled bending vibrations are ob- 
tained to reasonable accuracy, numerical data used 
in references [4, 6] are employed in the second 
order finite-difference method, and the natural 
frequencies and mode shapes are obtained. Although 
not presented here, the agreement between the re- 
sults obtained in this manner and those presented 
in reference [6] is found to be extremely close. 

Pretwisted Tapered Blades 

Among the recent contributions to pretwisted, 
tapered blade vibration analysis, first order 
f initp-differences were used by Carnegie et.al. 

[19, 20]; the Galerkin method was used by Rao [21]; 
finite element method was used by Gupta and Kao 
[22]; and the Reissner method was employed by 
Subrahmanyam and Rao [6]. Numerical data used in 
references [19, 20] are employed in the present 
Investigation, and several cases of tapered and 
pretwisted blades are solved. A typi al set of 
results, for a 45* pretwisted blade with a breadth 
taper of -0.25 and a depth taper of 0.75, are pre- 
sented in Figure 2 in order to represent the rela- 
tive convergence trends of first and second order 
theories. It has been observed [14] that, the 
finite-difference method clearly yields a lower 
bound for the uncoupled vibration cases, but for 
coupled bending-bending vibrations, the convergence 
is seen to be from above for certain coupled modes. 
As can be seen in Figure 2, the first order theory 
has oscillatory convergence for the first coupled 
mode while the second ordpr theory produced an 
upper bound. However, the higher modes are clear 
lower bound solutions. Extensive convergence 
studies were made for a variety of blades having 
various breadth-to-depth ratios and pretwist angle 
combinations discussed in [6, 19 to 21] by using 
both the first and second order finite-difference 
approaches. It appeared that for rectangular cross 
section blades, 'hen the vibratory mode was such 
that the maximum component deflection was in the 
flexible plane, the corresponding coupled frequency 
converged from above and when the maximum component 
deflection was in the stiffer plane, the associated 
coupled frequency yielded a lower bound. For 
blades with a square cross section at the root, 
the flexural rigidities are equal in both the 
planes and pretwist does not bring in any coupling. 
Coupled vibration can only occur if such a blade 
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has unequal breadth and depth tapers and also 
possesses a non-zero pretwist. In such cases, the 
convergence yields upper or lower bound depending 
upon the unequal flexural rigidities due to taper, 
the degree of coupling, and the mode of vibration. 
In general, for pretwisted tapered beams having 
square cross section at the root presented here, 
the convergent has been observed to be lower bound 
except for the first coupled mode. 

The authors have extended this work more 
recently to analyze rotating blades of asymmetric 
cross section, without pretwist, but which perform 
coupled vibrations having coupling between bending 
in one plane and torsion. Those results, although 
not shown here. Indicated that the second order 
finite difference method predicts close lower bound 
solutions for the coupled bending-torsion modes. 

8y comparing the trends observed for the un- 
coupled vibration cases, coupled bending-bending 
vlbravort' of pretwisted blades and the coupled 
bending-t orsion vibrations of untwisted asymmetric 
cross section blades under rotation, it is con- 
cluded <t .‘. t the coupling between the bending In 
two planes, brought forth by pretwist, is respon- 
sible for the occurrence of both upper bound and 
lower bound convergence in the finite-difference 
approaches depending on the mode of vibration. 

Figures 3 to 5 show some typical results pro- 
duced by the second order central difference method 
with n m 30. For certain values of breadth and 
depth taper, the effect of varying pretwist is 
shown in Figure 3. Figure 4 shows the effects of 
depth taper for a given value of breadth taper and 
pretwist. Figure 5 shows the coupled bending- 
bending mode shapes of a 90* pretwisted blade hav- 
ing a breadth taper of -0.25 and various values of 
depth taper. Further results, obtained for 
B— 0.5 ; e«0.5 and various values of a and v, are 
not shown here, but all these results are in good 
agreement with the theoretical and experimental 
values available in the literature [19. 20]. The 
effects of orptwist, breadth taper and depth taper 
in producing coupling of the modes are discussed 
in detail in references [6, 10, 17, 19]. 

Relative Computational Efficiency 

In order to evaluate the relative computa- 
tional efficiency of the first and second order 
finite difference methods, the average CPU time 
required by each method was determined for n « 30. 
The methods require nearly the same amount of CPU 
time; the maximum variation between the two is 
about +1.7 percent. 

CONCLUDING REMARKS 

The second order finite-difference method has 
been applied to determine the coupled bending- 
bending frequencies and the mode shapes of pre- 
twisted and tapered cantilever beams. Simple 
recursive relations have been used to eliminate 
the fictitious stations outside the beam domain by 
making logical extensions from the first order 
theory. The present approach is shown to produce 
accurate natural frequencies and mode shapes. In 
the course of this study, several conclusions have 
emerged. 

1. For the same mesh si 2 e (step length h), 
the second order finite difference method produces 
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natural frequencies with greater accuracy than the 
first order theory. The convergence of the lower 
mode frequencies is much more rapid for second 
order central differences than for first. 

2. Second order theory yields natural fre- 
quencies and mode shapes to accuracies of practical 
Interest with relatively coarse mesh. Further, 
extrapolation procedures are not necessary to ob- 
tain accurate results. 

3. Results presented here Indicate that the 
second order theory produces close lower bound 
solutions for uncoupled mode of vibration. For 
the modes having coupling between bending In two 
planes, either close lower or close upper bound 
solutions are obtained. The probable conditions 
under which an upper bound solution Is obtained 
varied with the breadth to uepth ratio, pretwist, 
taper ratios and the extent of coupling between 
the two component modes. 

Majority of short comings encountered In the 
earlier investigations such as the necessities of 
using integer and half integer stations, of trans- 
forming the equations of motion into a set of first 
order equations or of using backward or forward 
differences at the boundaries to avoid fictitious 
stations, etc., are eliminated by the present im- 
proved theory. Extension of the method presented 
here using plate theories may prove beneficial, 
since the strong convergence characteristics would 
reduce computational time to a considerable extent 
in two dimensional cases. 
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TABLE I. - CONVERGENCE PATTERN OF FLEXURAL FREQUENCIES OF UNIFORM CANTILEVER 
BEAM USING FIRST ORDER AND SECOND ORDER CENTRAL DIFFERENCE SCHEMES: 



NONDIMENSIONAL FREQUENCY, pi 


n 

I Mode 

11 Mode 

III 

Mode 


1st order 

2nd order 

1st order 

2nd order 

1st order 

2nd order 

: 

3.4021 


1° .870 

21.204 

45.334 

54.421 

6 

3.4359 


19.107 

21.553 

49.370 

57.320 

10 

3.4866 

3.5148 

21.134 

21.934 


60.774 

11 

3.4916 

3.5162 

21.284 

21.960 

57.423 

6) .014 

IP 

3.4955 

3.5153 

21.340 

21.977 

58.063 

6. .179 

15 

3.5029 

3.5158 

21.623 

22.006 

59.318 

61.443 

1? 

3.5057 

3.5158 

21.713 

22.015 

59.827 

61.526 

18 

3.5069 

3.5158 

21.750 

22.018 

60.024 

61.555 

20 

3.5086 

3.5159 

21.801 

22.023 

60.334 

61.595 

23 

3.5104 

3.5159 

21.857 

22.027 

60.660 

61.632 

24 

3.5108 

3.5159 

21.872 

22.028 

60.744 

61.640 

<5 

3.5113 

3.5159 

21.884 

22.029 

60.817 

61.64 7 

2S 

3.5125 

3.5160 

21.923 

72.031 

61.041 

61.666 

30 

3.5127 

3.5160 

21.930 

22.031 

61.083 

61.669 

Exact 

r 

3 .5160 

22. 

0345 

61 .6973 

soluti on 







Percent 







error 







based on 








-0.094 

0.000 

-0.474 

-0.016 

-0.996 

-0.046 


TABLE II. - COMPARISON OF FREQUENCY PARAMETER RATIOS, y- , OF PRETWISTED UNIFORM BLADES 


V c o 

ratio 

— 

V 

Mode 

not.^er 

First order finite- 
difference method: 
extrapolated from 
n - 10. 15. 20 [19] 

Matrix displacement 
method [18] extra- 
polated from 

Second order 
f i ni te-dif ference 
method: n - 30 

without extrapolation 

3 and 4 
elements 

8, 10 and 16 
el ements 

2 

30* 

i 

1.0059 

1.0049 

1 .0015 

1.0053 



2 

3.9166 

3 9168 

3.9132 

3.9139 



3 

40.3667 

40.3834 

40.5648 

40.3745 



4 

148.8397 

149.0163 

149.3350 

148.6643 

| 

90' 

i 

1 .0425 

1 .0429 

1 .0444 

1.0434 



2 

3.4020 

3.4037 

3.3935 

3.3968 



3 

48.6700 

48.6587 

48.65H7 

40.7710 



4 

113.4711 

113.9809 

113.3201 

113.0779 

4 

30“ 

i 

1 .0044 

1 .0060 


1 .0081 



2 

13.9468 

13.9524 

— 

13.8957 



3 

45.3107 

45.3285 



45.4342 



4 

278.3168 

280.9137 

— 

277.6319 



5 

693.4500 

— 

— - 

693.9047 




















NATURAL FREQUENCY, Hz 



Figure Z - Convergence pattern of coupled flexwial fre- 
quencies. (b a /c 0 * 1, (5 * -0. 25, 6 ■ 0. 75, r ’ 45°). 
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Figure 3. - Frequency parameter 
ratio against pretwist angle for 
a blade with p • -a 25, 6 • 0. 75 
and tyc 0 -l. 
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